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Abstract 

In this paper, we derive coupled envelope equations modeling the growth of stimulated Raman 
scattering (SRS) in a multi-dimensional geometry, and accounting for nonlinear kinetic effects. 
In particular, our envelope equations allow for the nonlinear reduction of the Landau damping 
rate, whose decrease with the plasma wave amplitude depends on the rate of side-loss. Account 
is also made of the variations in the extent of the plasma wave packet entailed by the collisionless 
dissipation due to trapping. The dephasing between the electron plasma wave (EPW) and the 
laser drive, as well as the self-focussing of the plasma wave, both induced by the EPW nonlinear 
frequency shift, are also included in our envelope equations. These equations are solved in a 
multi-dimensional geometry using our code dubbed BRAMA, whose predictions regarding the 
evolution of Raman reflectivity as a function of the laser intensity are compared against previously 
published PIC results, thus illustrating the ability of BRAMA simulations to provide the correct 
laser threshold intensity for SRS, as well as the right order of magnitude of Raman reflectivity 
above threshold. 
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I. INTRODUCTION 



The present work is part of the effort that has been undertaken for about one decade by 
several groups in order to understand, and efficiently model, the impact of nonlinear kinetic 
effects on the growth of stimulated Raman scattering (SRS). All the recent modelings we 
know of, including the one presented here, rely on the hypothesis that the wave electric 
fields may be written in terms of a slowly varying envelope and a rapid phase. Using 
this approximation, Vu et al. developed a reduced Particle In Cell (PIC) code, described 
in Ref. [l| (and a further simplified PIC-like scheme was even derived by Hur et al. in 
Ref. j^), which they used in Ref. {3] to show that Raman reflectivity could reach much higher 
levels than linear theory would predict. This so-called "kinetic inflation" was confirmed 
experimentally by Montgomery et al. in Ref. as well as by full kinetic simulations (see 
for example Refs. , and was attributed to the nonlinear reduction of the Landau 

damping rate of the electron plasma wave (EPW), following the original idea of O'Neil 7]. 
Such results triggered several theoretical developments, for example those of Refs. js-12], 
revisiting O'Neil's work and aiming at a formulation of the collisionless damping rate general 
enough to apply to a realistic physics situation, such as SRS. 

Another important nonlinear kinetic effect, widely discussed this last decade, is the non- 
linear frequency shift of the EPW, which may induce a dephasing between the laser drive 
and the plasma wave and, therefore, hamper the growth of SRS, as argued for example 
in Ref. p]. Moreover, in more than one dimension (1-D), the frequency shift usually has a 
transverse profile making the EPW phase velocity smaller close to the wave axis of propaga- 
tion than away from it, causing the so-called "wave front bowing" evidenced numerically in 
Ref. jl3j. In practice, this entails a transverse component in the plasma wave number that 
tends to make the EPW self-focus. Just like for the nonlinear collisionless damping rate, 
several recent theoretical papers (for example Refs. Is, 14, 15]) revisited the original works 
on the nonlinear frequency shift of an EPW, Refs. la, Il7 |. 

The aforementioned theoretical works were mainly motivated by the will to develop en- 
velope codes running fast enough to address such a large scale system as a fusion hohlraum 



18j, and accurately accounting for nonlinear kinetic effects. Such envelope codes were de- 
scribed in Refs. 19, 20] and their predictions were compared with those of 1-D kinetic 
simulations. In this paper, we go further in that direction, and show what we believe are 
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unprecedented comparisons between results from PIC simulations and from our envelope 
code, BRAMA, for Raman reflectivity in a two-dimensional (2-D) and a three-dimensional 
(3-D) geometry. 

In BRAMA, the modeling of SRS rests on the hypotheses that the plasma is collisionless 
and non relativistic, and that the waves are purely sinusoidal. The non-relativistic approx- 
imation is clearly valid for the laser intensities used in our simulations. Moreover, for the 
range of parameters considered in this paper, the main effect of collisions, which is to par- 
tially restore Landau damping (see Ref. [l9[ for details), is overcome by that due to the 
transverse electron detrapping. The sinusoidal assumption is vindicated a posteriori by our 
results showing that SRS saturates at EPW amplitudes small enough for the effects of an- 



harmonicity to be negligible. This allows us to use in BRAMA the values derived in Ref. [14 1 
for the frequency shift, 5u p , which are in very good agreement with those inferred from 1-D 
Vlasov simulations of SRS, and which rest of the hypothesis that the electron motion is 
nearly adiabatic. The latter approximation is valid when l± > A;, where l± is the transverse 
extent of the plasma wave and A; the laser wavelength (see Ref. [ll| and Appendix EJ, a 
condition we were careful to meet in our BRAMA simulations. The dephasing between the 
plasma wave and the laser drive induced by 5u p is accounted for within the hypothesis of 
self-optimization, which is recalled in Section ITTTI . and which also is supported by 1-D Vlasov 
results (see Ref. j^lj]). As for the self-focussing induced by 5oo p , it is explicitly contained in 
our envelope equations. 

The nonlinear expression we use for the EPW collisionless damping rate follows from 
the theoretical developments of Ref. [l^ . and allows for transverse electron trapping and 
detrapping. This value for the damping rate was checked in Ref. 12J to agree with results 



obtained from 1-D Vlasov simulations of SRS, and was shown in Ref. [23j to be consistent 
with that of Yampolsky and Fisch [l0| over a finite range of wave amplitudes. In BRAMA, 
we also account for the collisionless dissipation of the plasma wave induced by trapping, 
which entails no damping but alters the extent of the plasma wave packet (both in the 
longitudinal and transverse directions) due to a nonlinear and non local change in the EPW 
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29[ to be in very good 



group velocity, whose theoretical estimate was checked in Refs. 
agreement with results from 1-D Vlasov simulations of SRS. 

We are therefore using in this paper a nonlinear kinetic modeling of SRS which we checked 
very carefully using 1-D Vlasov simulations. However, other nonlinear kinetic effects than 
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those previously cited exist, for example the growth of sidebands evidenced in Refs 



321 ]. and may have important implications in SRS. These are nevertheless not accounted for 



in BRAMA, and, in spite of the imperfections of our modeling, we find it interesting to 
discuss its ability to make relevant predictions as regards Raman reflectivity. This is the 
purpose of the present paper to do so by showing comparisons between BRAMA results and 
the PIC simulation results obtained by Yin et al. in Refs. jf], Q] . 

This article is organized as follows. In the next Section, we derive the envelope equations 
solved in BRAMA, and discuss them physically. Section II I II is devoted to the comparisons 
we made between the results from PIC and BRAMA simulations regarding the evolution of 
Raman reflectivity as a function of the laser intensity. Section HVl summarizes and concludes 
this paper. Moreover, supplementary material and discussions are given in the Appendices. 



II. WAVE EQUATIONS 

In this paper, Raman scattering is studied within the context of the so-called three-wave 
model, wherein the total electric field reads 

Em = -z^le^+MO _ ^| e ^ in + fc| c *S» + c . c ., (i) 

where (p l ™ = k l ^x — ui l ™t, with w — p,l,s respectively for the plasma, laser and scattered 
wave, and (pf n = ip l p in + ip 1 ™. Moreover, (cuj)°) 2 = oo pe + (fcj^c) 2 , c being the speed of light in 
vacuum and u pe the plasma frequency. cuf n is the vacuum laser frequency, w. 
are such that they maximize the linear SRS growth rate derived in Ref. 27 
accounts for the effect of the frequency shift, 5u p = —d t 5ip p , while the £ w 's are slowly varying 
envelopes, {S^dxSwl <C |fc^ n | and {S^dtSwl "C |w^ n |. Moreover, as shown in Appendix [Aj 
the polarization x p is very close to the direction of the wave number, k p = k l p m x + V5ip p , 
because the plasma wave is nearly electrostatic. As for the laser and scattered waves, we 
account for the fact that they are not exactly polarized along the y direction, and derive in 
Appendix [A] the relation between their x and y components. 



rile fc; m and u p in 
I. As for 5tp n , it 



A. Equation for the plasma wave 



In order to derive the envelope equation for the plasma wave amplitude, we use the 
same procedure as in Ref. 22J which consists in first using the results from Whitham's 
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variational approach [28[ , and in then adding "by hand" the effects of dissipation and of the 
electromagnetic drive, which are not accounted for in Whitham's theory. This procedure 
can be justified a posteriori by going through Vlasov-Gauss equations as was done in the 



Appendix of Ref. 



29), or in Ref. Hi 



This program will however not be pursued here, and 
we prefer refer the reader to Ref. [22j where our procedure appears quite naturally. 
Hence, we first start by using Whitham's variational approach, which yields 

3 2 L 



dtdco. 



- V. 



0. 



(2) 



where u p = u p m —dtS(f p is the nonlinear frequency of the plasma wave. We henceforth assume 
that £ p is real so that the Lagrangian density is L = J^ p ds'Ld£' } where ds/L = is the 
dispersion relation of a freely propagating EPW, at 0-order in the variations of £ p . At this 
order, the plasma wave may clearly be considered electrostatic, so that d^L = (1 + Xr)£', 
where \ r is the adiabatic approximation of the real part of the electron susceptibility. It 
assumes the same values as in 1-D, derived in Ref. 11] . and only depends on the plasma 
wave number through its modulus, k p , since the EPW may be considered electrostatic and 
the electron distribution function is isotropic. Then, using L = / p [1 + Xr£'} dS' and the 
consistency relation d t k p = — Vw p , Eq. fl2]) yields, 



dxr d£ P 
duip dt 



d^CrdSp Sp 



dt dul p dk n pA v dkl 



0. 



(3) 



where x p = k p /k p and, therefore, d Xp £ p ee {k p .V£ p )/k p . Due to the consistency relation 
dtk p = —Vu p , a transverse profile of the EPW frequency entails a transverse component in 
k p , usually directed towards the wave axis of propagation. Hence, in Eq. (j3]), we do account 
for the self-focussing induced by wave front bowing, as discussed in Ref. 13]. Note, however, 
that making use of an adiabatic dispersion relation does not allow for the phase modulation 
induced by the transverse dependence of the wave amplitude. Hence, we do not account for 
diffraction-like effects in Eq. (j3J), and this point is discussed in detail in Appendix iBl 

Now, since the nonlinear frequency shift only depends on £ p , dtu p is directly proportional 
to d t £ p , and, actually, £ p d^ p Xr9t^ P *C d Up x r dt£p because 5u p u p . Moreover, within the 



paraxial approximation, k^^ k px , which yields V.x p ~ k p 1 V±.k p 



k p 1 + d z k 



Finally, as shown in Ref. [2jJ, wherever SRS is effective, 5k p w 5uj p /v gS) where v gs ~ c is the 



group velocity of the scattered wave, which makes Xp.Vk p directly proportional to d Xp £ p . 
Then, because 5u p u p and the EPW phase velocity, t>^, is much less than v gs , the term 



proportional to x v SJk p in Eq. 03) is much smaller than that proportional to d Xp £ p . Conse- 
quently, since in this paper we are mainly interested in the kinetic modeling of stimulated 
Raman scattering, we simplify Eq. (j3J) in 



dxr dS p dxr 



duj p dt 



We now need to add to Eq 
discussed in Ref. 



dk r , 



BE £ - - ' 



0. 



(4) 



the terms accounting for collisionless dissipation. As 



221 ]. two different effects need to be allowed for, a Landau-like damping 



and the irreversible exchange of energy between the wave and the electrons induced by trap- 
ping. The expression we use for the collisionless damping rate follows from the theoretical 
developments of Ref. 12J, which would yield, if all electrons had the same transverse velocity 
v±, 

p(v±) = v Xm u [y 3D (v ± )} , (5) 

where z/ lin in the (linear) Landau damping rate [30[, Y 3D = J"* w B (x// — v^t', x± — v±t', t')dt', 
ub = a/ eSpm/kp being the bounce frequency, and %(Xsd) is a Heaviside-like function, 
U(Y 3D ) « if F 3 D > 6 and U(Y 3D ) w 1 if Y 3D < 6. 

As for the rate of dissipation induced by trapping, it is clearly proportional to the rate 
of electron trapping, which, if we neglect the small change in the EPW phase velocity, is 
just proportional to the wave growth rate, calculated in the wave frame, and as seen by 
the electrons. For electrons with transverse velocity vj_, the rate of dissipation induced by 
trapping is therefore proportional to d t £ p + v<f,d Xp £ v + v±.V±£ p . Quite naturally, and as 



shown in Refs. 
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22], in the envelope equation for the EPW the prefactor of the latter 



expression is —d Up ^, which is the contribution of the trapped electrons to —d Up Xr ■ More 
precisely, xjL r (r, £) is the contribution to the real part of the electron susceptibility of 
those electrons with transverse velocity vj_ which have been trapped by the EPW at time 
t' < t, and have completed at least one trapped orbit (i.e., Y 3 b> ^ 6). Hence, x* r is a non 
local function of the plasma wave amplitude, which mainly depends on the maximum EPW 
amplitude experienced by the electrons with transverse velocity v± (see Ref. 22| for more 
details). 

We therefore conclude that, when all electrons have the same transverse velocity, Vj_, 
allowing for collsionless dissipation leads to the following envelope equation for the plasma 
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wave amplitude 



env 
r 



dUJr 



dt 



dXr 

dk 



' dujp 
dco P X 



where we have introduced d Wp X 



v 

env 



Refs. 
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d£ 



V 



dx 



tr 



(hp dk p 2k p ' ^"" p dujp^ ^ p ^' 
du> P x) r , whose expression may be found in 



20] . When the transverse electron motion is not infinitely cold, one needs to 



average Eq. ([6]) over the distribution, f(v±), of transverse velocities, which we approximate 
by the unperturbed one (assumed to be a Maxwellian) , since the quiver transverse velocity 
induced by the laser field is much less than the thermal speed for the parameters investigated 
in this paper. Using the facts that d^ p Xr does not depend on v±, that f(v±) is isotropic, and 
that d kp Xr = {-2/k p )xr ~ v+d^Xr ~ (2/kp) - v+d^Xr (see Ref. fl for details), averaging 
Eq. dSJ) over v± yields 



env 
3D 



du Ti 



d£ - £ p - -> 
— f + v g V£p + v Q -f-V±.kp + u 3D £ p 
dt 2k p 



0. 



with 



- y* 



J 9y,z 



V 



"3D 



I'd, 



Vcj> 



V<f, 



h f) v env 
2 

h f) v env 
9u) v Xr 



y,Z I d Wp X? V f(v±)Vy,zdVj 



k 



d. P Xt% 



X v env u a v env ' 

I f(v±Hv±)d UpX env ^±)dv± 



(7) 

(8) 

(9) 
(10) 

(11) 



dupXtZ 

where d Up xtj) = / f(v±)d LJp xl nv dv±. Note that, in order to derive v g , we made use of the 
approximation, v± « v y y + v z z. Note also that Xm an d du P Xl nv are non local functions of 
the EPW amplitude, which entails the dissipation of the electrostatic energy. More precisely, 
the expression we find for v gx , and for the second term of Eq. ([9]) for v gy and v gz , respectively 
lead to the longitudinal and transverse shrinking of the plasma wave packet in the strongly 
nonlinear regime, while maintaining the EPW amplitude constant along the characteristics, 
which automatically reduces the total electrostatic energy (see Ref. 22] for details). As 
for the first term in Eq. OH]), it allows for the self- focussing induced by wave front bowing. 
Moreover, it is quite clear from Eq. ( TTTj) that we account for side-loss in our expression for 
the collisionless damping rate v 3 r>, and, from Eq. flH]), that the latter would decrease less 
rapidly as a function of the EPW amplitude for a larger rate of side-loss (narrower wave 
packet or larger electron thermal velocity). 



In order to complete our envelope equation for 8 P we now need to account for the effect 
of the ponderomotive drive. From Eq. (IA.8j) of Appendix |XJ it is clear that, up to terms 
at second order in k y /k, which we neglect, the ponderomotive force is the same as in 1-D. 
The right-hand side of the envelope equation for the EPW, accounting for the laser drive, 
is therefore exactly the same as the one derived in Ref. [^(J, so that, if we denote 

E p = S p /2, (12) 
E s = (£ s /2)e i5 ^y.y s , (13) 
E l = (£,/2)y.y,, (14) 

Tp = muj^^Xm' (15) 

we find 

dE E 

+ v g VE p + v Q ^V ± .k p + u 3D E p = T p EiE* s . (16) 

Note that the solution of Eq. f JTB"]) is not necessarily real as E p should be. Actually, enforcing 
the reality of E p would give rise to technical difficulties when numerically solving the EPW 
envelope equation, and would actually only make sense if we were able to perfectly calculate 
5ip p , which is not the case. We therefore relax the constraint that E p needs to be real, 



and, although our procedure is not totally rigorous, it proved in Refs. [20|, [2lj to provide 
results in very good agreement with those of 1-D Vlasov simulations. This is because the 
only phase that has an impact of the growth of SRS is not that of the EPW itself, but the 
dephasing between the plasma wave and the laser drive, which we estimate by making use 



of the self-optimization ansatz detailed in Ref. 21| . and which is recalled in Section 



B. Equations for the laser and scattered waves 

The equations for the laser and scattered waves are derived from Maxwell laws by using 
the paraxial relation derived in Appendix El E x « (i/k)d y E y , between the components 
of the electric field. More precisely, plugging Eq. ( 1A.8I) into Eq. ( 1A.4I) and making the 



usual approximations which consist in neglecting E z , the space variations of k, the second 
derivatives in x and t as well as space derivatives of order larger than 3, yields 

[2ikd x + 2i{u/c 2 )d t + (u 2 /c 2 ) -k 2 + A J E y = fi (d t - iu)j y , (17) 
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where Aj_ = d 2 + d 2 . As is well known, fio(d t — icu)j y is the sum of the linear response to 
E y , which yields u 2 e E y , and of the "Raman" term proportional to the density fluctuation, 
8 ne i(<P P n +s<p P ) + CC-) induced by the EPW (see Ref. |20j for details). Since the relation 
ik p Sp ~ —5ne/Eo still holds, the "Raman" term is the same as in 1-D, so that one recovers 
the equations derived in Ref. [20|, except for the additional term proportional to A^E y . 
Hence, the envelope equation for the laser wave amplitude, Ei, is 



[d t + v gl d x - i(v g i/2ki)A ± ]Ei = TiEpEg 



(18) 



where v g i = kic 2 /uji and Ti = ek p /2mu s . As regards the envelope equation for the scattered 
wave amplitude, E s , we need to account for the extra phase 5(p p in its definition, which 
yields, 



d t + v gs d Xp + i(5u p - v gs 5k p ) - i(v gs /2k s ) (A ± + iV±.k t 



where v gs = k s c 2 /u s and T s = ek p /2mu!i. 



III. COMPARISONS WITH PIC SIMULATION RESULTS 



e s = r s EiE;, (19) 



by our code, BRAMA, 



201 ]. with the additional 



The previously derived envelope equations are solved numerically 
using the standard operator-splitting approach detailed in Ref. 
hypothesis of self-optimization which consists in making the dephasing term (5u p — v gs 5k p ) 
in Eq. ( fl9l) as small as possible. More precisely, we calculate 5k p from the consistency 
relation d t Sk p = —V5ujp, and, when we find that Sk p thus calculated may be larger than 
Sujp/v gs , we just fix it at 5<jj p /v gs . This procedure proved to provide results in very good 
agreement with those of 1-D Vlasov simulations of SRS in Ref. [20 ] . 

Using BRAMA we aim at deriving the laser threshold intensity for SRS, as well as an 
order of magnitude for Raman reflectivity above threshold. In order to test the ability of 
our envelope code to do so, we compared its predictions against the 2-D and 3-D simulation 
results published by Yin et al. in Refs. js, Q], some of them we double-checked using our 
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own PIC code CALDER _ 

Let us start with the 2-D simulation results of Ref. js], corresponding to a plasma with 
electron temperature T e = 700 eV and electron density n e /n cr = 0.036. From the left end of 
this plasma is injected a laser with wavelength = 0.527 /im, which propagates along the x 
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FIG. 1: (Color online) Raman reflectivity, R, for a 2-D geometry, as obtained in Ref. [6fl using PIC 
simulations (red circles), and derived from BRAMA simulations with rj = 1CP 5 (green squares) or 
r) = 10~ 10 (blue diamonds). 

direction and whose envelope at best focus, centered in simulation region, is a 2-D Gaussian 
whose intensity varies as exp(— z 2 /w 2 ) with a waist w = 2.58 /im. In BRAMA simulations 
is also launched, from the right end of the plasma, a seed laser which is not focussed, whose 
peak intensity, I s , is related to that of the laser by I s = r]Ii aser , and whose wavelength is 
X s ~ 0.684 /im (which maximizes the linear SRS growth rate for the chosen temperature 
and density). Then, the plasma wave generated by Raman is such that Ai m A£> ~ 0.34, 
where Ad is the Debye length. The simulation domain has size 100 /mix 14 /im, which, in 
BRAMA simulations, we discretize using 400 points in the x direction and 64 points in the 
y direction. With BRAMA, it takes about 10 minutes on one processor to simulate 1 ps 
of laser-plasma interaction. As may be seen in Fig. (TJ when rj = 10~ 5 , the same trend is 
observed as regards the intensity dependence of Raman reflectivity calculated using either 
PIC or BRAMA simulations, which both predict that Raman scattering should not be very 
effective when 7i ase r ^4x 10 15 W/cm 2 (although the onset of SRS is more marked in PIC 
than in BRAMA simulations). 
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FIG. 2: (Color online) Raman reflectivity, R, for a 3-D geometry, as obtained in Rcf. 26] using 
PIC simulations (red circles), and from BRAMA simulations with rj = 10 -5 (green squares). 



Now, one may wonder whether the good agreement found for the SRS threshold is not 
just fortuitous, due to our particular choice for rj, since in BRAMA Raman scattering results 
from the optical mixing of two counterpropagating lasers, while in PIC simulations it starts 
from the numeric noise. In order to address this issue, we strongly decreased rj (since we 
cannot greatly increase it without missing the linear regime), down to rj — 10~ 10 , which 
is way below the level of electromagnetic noise measured in our own PIC simulations with 
CALDER. As shown in Fig. [I], with such a low value of i] the threshold intensity for Raman 
scattering is Ith ~ 7 x 10 15 W/cm 2 . Hence, the dependence of with respect to 7] is very 
weak, a decrease of r\ by 5 orders of magnitude leads to an increase in I h by a factor less 
than 2! This shows the little impact of the choice of rj on our results, and therefore confirms 
our ability to predict the SRS threshold using BRAMA. 

As regards Raman reflectivity above threshold, BRAMA results systematically overesti- 
mate the PIC ones, especially when /i aS er = 10 16 W/cm 2 . For this particular intensity, we 
find that, as SRS reflectivity increases, the EPW experiences a strong self- focussing due to 
wave-front bowing, which should limit Raman scattering. More precisely, when R < 1%, 
self-focussing is not yet effective, so that our modeling relying on the paraxial approxima- 
tion should be valid, and we can therefore predict that, when Ji aser = 10 16 W/cm 2 , Raman 
reflectivity should reach at least 1%. When R « 6%, the EPW is strongly focussed, its 
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the transverse extent is significantly less than the laser wavelength, and we have to stop 
the simulation. Indeed, such a result is clearly not physical and, although we are currently 
not able to prove it, it is most likely that Raman should saturate by some means (usually 
due to the growth of sidebands) before R reaches 6%. We therefore conclude that, when 
-fiaser = 10 16 W/chi 2 , SRS should saturate at a value of R between 1% and 6% (and it is 
noteworthy that, in the PIC simulations of Ref [6], the first maximum of SRS reflectivity 
is at R « 6% when /i aS er = 10 16 W/cm 2 ). Hence, for the time being, BRAMA simulations 
cannot provide an extremely accurate estimate for the long time averaged Raman reflec- 
tivity, but can nevertheless predict what should be its order of magnitude, typically a few 
percents when 7i aser = 10 16 W/cm 2 , which is consistent the PIC results of Ref. 6| yielding 
(R) ~ 1.5%. It is also noteworthy that, as shown in Fig. [U our results for R above threshold 
do not depend on our choice for rj. 

The same conclusions as before hold for the 3-D simulation results plotted in Fig. [2j 
which correspond to a plasma with electron temperature T e = 4 keV and electron density 
n e /n cr = 0.14. The laser wavelength is A/ = 0.351 /im, its waist is w = 1.4 /im, and it 
is focussed in the center of the simulation box whose size is 35 x 6 x 6 /im 3 . In BRAMA 
simulations, the seed wavelength is A s ~ 0.626 /im (the corresponding EPW is such that 
kp n \o ~ 0.32). The simulation domain is discretized using 100 points in the x direction and 
32 points in the y direction. With BRAMA, it takes about 1 minute on 128 processors to 
simulate 1 ps of laser-plasma interaction. As may be seen in Fig. |2} and as for the previous 
2-D simulations, we observe the same trend regarding the evolution of R as a function of 
/laser using PIC or BRAMA simulations. In particular, Raman reflectivity is expected to 
be very small when Ji ase r < 10 16 W/cm 2 , and to be of the order of a few percents when 
-fiaser = 5 x 10 16 W/cm 2 , although, for this particular intensity, the BRAMA prediction for 
R overestimates that of the PIC simulation. 

Using BRAMA, we actually reran all the simulations presented in Refs. I^Q, and always 
found results similar to those shown in Figs. [Hand 12 We may therefore conclude that, both 
in 2-D and in 3-D, BRAMA is able to correctly estimate the SRS threshold and that, above 
threshold, it provides the correct order of magnitude for R. Therefore, using BRAMA one 
should be able to conclude about the effectiveness of Raman scattering, at least for a single 
laser speckle. 
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IV. CONCLUSION 



In this paper, we derived coupled envelope equations modeling stimulated Raman scat- 
tering in a multidimensional geometry, and allowing for such nonlinear kinetic effects as the 
frequency and wave number shift of the electron plasma wave, the nonlinear reduction of its 
collisionless damping rate, and the collisionless dissipation induced by trapping. 

Our theoretical predictions for the aforementioned nonlinear kinetic effects were shown 
to be Jiii ver y g ood agreement with results from 1-D Vlasov simulations of SRS in 



Refs. |12l . Il4j . |21| . We are therefore using here a nonlinear kinetic modeling of SRS that 
we previously checked very carefully, and that we merely generalized to allow for a multi- 
dimensional geometry so as to derive envelope equations which we solved using our envelope 
code BRAMA. In order to test BRAMA predictions as regards Raman reflectivity, we took 
advantage of the numerous 2-D and 3-D PIC results published by Yin et al. in Refs. jf], Q]. 
When making comparisons which such PIC results, we had to be careful that, in BRAMA, 
we simulate the optical mixing of two counterpropagating lasers, while SRS starts from nu- 
merical noise in PIC simulations. In particular, we checked that the good agreement we 
found between the PIC and BRAMA predictions as regards the laser threshold intensity, 
Ith, was not fortuitous, not due to the particular choice we made for the seed intensity, I s , 
by showing the very weak dependence of I t h with respect to I s . 

Clearly, in BRAMA, we do not account for all possible nonlinear kinetic effects and, in 
particular, not for the growth of sidebands, although this is something we are currently 
working on, following the lines of Ref. j^J. Consequently, we cannot make any claim re- 
garding our ability to predict the long time averages of Raman reflectivity. We nevertheless 
interpret the strong EPW self-focussing we observe in BRAMA simulations as a natural limit 
for Raman growth, which allows us to predict an order of magnitude for SRS reflectivity 
above threshold, which appears to be consistent with that of the PIC simulations. Hence, 
although it is clear that the present paper is certainly not the last word on the nonlinear 
kinetic modeling of SRS, we may nevertheless conclude that, in its present state, BRAMA 
appears to be a very powerful tool to predict Raman effectiveness, at least for a single laser 
speckle. 
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Appendix A: Discussion on the nearly electrostatic, or electromagnetic, nature of 
the waves. 



Let us consider an electric field, £ , which we write 

£ = [E x x + E y y + E z z\ e^ kx -^ + c.c. (A.l) 
From Maxwell equations, we know that 

V 2 £- c- 2 d 2 t £ = fi dj + V[V.£], (A.2) 
where j is the current density. Projected on the x, y, and z directions, Eq. flA.2j) yields 

[d 2 + d 2 + c~ 2 (lu 2 - d 2 + 2iud t )] E x = fM (d t - iu)j x + (d 2 xy + ikd y )E y 

+(d 2 xz + ikd z )E z (A.3) 
[d 2 x + d\ y + 2ikd x -k 2 + c' 2 (to 2 - d 2 + 2iud t )] E VtZ = n Q {d t - iu)j y , 2 + ikd VjZ E x 

~) 2 F _|_ r) 2 P 
'xy^x ' u yz Ill z,y 



+d 2 E x + d 2 E z , y . (A.4) 



1. Nearly electrostatic nature of the plasma wave 

The nearly electrostatic nature of the plasma wave is best viewed by writing Eq. (1A.2|) 
in Fourier space, which yields, along the y or z direction, 

[W 2 /C 2 - K 2 ] £y )Z + i[lQWjy )Z = -Ky )Z (K X £ X + Ky£y + k z £ . (A. 5) 

For a freely propagating wave, inowj y>z = —{w 2 /c 2 )£ yyZ) while for a driven wave one needs 
to account in ifi wj y for the additional term —(w 2 /c 2 )£d y , induced by the y component of 
the driving field. Now, because the electron thermal velocity is much less than the speed of 
light, w 2 /c 2 k 2 ~ k 2 . Moreover, from Eq. flA.7[) . £^ y ~ (n y / K x )£ dx , while \£d x \ "C \£ x \- 
Hence, we conclude that {w 2 / c 2 )\£d y \ is much less than \n y K x £ x \, and is thus negligible 
in Eq. flA.5j) . whose left-hand side may therefore be approximated by —K 2 £ y ^ z w —K 2 £ y , z . 
Moreover, since we only aim at modeling the nearly paraxial propagation of the plasma 
wave, which is essentially polarized along the x-direction, the right-hand side of Eq. flA.5[) 



is close to —K,y )Z K x £ x . Hence, 

£y,z ~ {l^y,z K x/ K )£x ~ {^"y,z/ K)£ x . (A-6) 

The significant Fourier modes of the EP W electric field are nearly electrostatic, and, there- 
fore, so is the total field. 
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2. Longitudinal component of the laser and scattered wave electric fields 



We now assume that E is essentially polarized along the y direction, which is the case 
for the laser and scattered waves. Then, at lowest order in the field variations, Eq. (|A.3[) is 



oj /c )E X = -ip uj x + ikd y E y . 



(A.7) 



Since —iujj x w u 2 e E x (if we neglect the density modulation), and provided that u and k are 



chosen so that uj 2 



pe 



+ k 2 c 2 , Eq. (TA~7j) yields 

E x « (i/k)d y E y . 



(A.8) 



As for E z , from Eq. flA.4j) we would find that it remains at very low levels, and we henceforth 
neglect it. 



Appendix B: Diffraction terms for the plasma wave 

Diffraction-like effects for the plasma wave are only addressed in the linear regime, where 
we slightly change our notations and write the EPW electric field -Eepw = E p e^ koX ~ UJot ^ +c.c, 
and the electron charge density p = p Qe l ( k o x -^ot) _|_ c c Then, if we introduce 

£oko(Epx + iE dx ) ' 

where, using the same notations as in Section[TTl E dx = ek EiE* / (moJiu s ) is the x-component 
of the driving field, Gauss law yields, 



(1 + X )E px - ^V.E P = -i X E dx . (B.2) 

Kq 

In order to calculate Xi we use a Fourier representation of the fields and write 

p = -ie J £{k, ^k.Ee^-^dkdu, (B.3) 

where E is the sum of the EPW and driving fields and, since each Fourier mode of the EPW 
field is nearly electrostatic, £ only depends on the wave number k through its modulus, k, 
as indicated in the integral of Eq. ( IB. 31) . This integral is approximated by making use of a 
Taylor expansion about k = k x and u = Uq, at first order in Sk x , 5u, 5k 2 and 5k 2 , 



£(k,u)k.E w £(k ,u) ) 



k E x + 5k. E 



k E x 



^(sk + 5 Sl+ + ^5u 

dk \ x 2k 2k J du 



, (B.4) 
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which yields 



i - f v.4 



_ , ae oe^ | , ae _^^ A±jE a (B 5) 

dk dx duo dt 2k dk px ' 



where we made use of the approximation E~ l d u E ps E~^d u E p for it = x,y,z,t. This approx- 
imation is vindicated in the regime of strong damping because, in this regime, E p and Ed 
are nearly proportional, while, when the damping is small, Ed <C E p and the approximation 
still holds. Plugging Eq. (|B.5|) into Eq. (1B.2j) . we find 

When the Landau damping rate, z/u n , is small compared to the plasma frequency, £(fco, w o) ~ 
£ r (fc , Wo) + «^ ^ r z/ lin , where £ r is calculated by making use of the adiabatic approximation. 
Hence, £ r assumes the same values as what we previously used for Xr m the variational 
approach of Section [III Thus, if we choose k and uo such that 1 + £ r (&o,^o) = 0, Eq. ( IB. 61) 
becomes 

[dt + v g0 d x - i(v g0 /2k )A ± + z/ lin ] E px = r p EiE*, (B.7) 

where v g0 = —dk ^ r /d Uo ^ r and where, as usual, we calculated the right-hand side of the 
envelope equation at lowest order in the variations of the field amplitudes, so that T p is 
given by Eq. (fl5|) . Is is quite clear that Eq. (1B.7j) is very similar to Eq. ([TBI) for the laser 
wave amplitude. Now, because v 9 q v g i )S , it is also quite clear that diffraction is much less 
effective for the plasma wave than for the laser or scattered wave. Since, as regards SRS, 
only the relative phase between the plasma wave amplitude and that of the ponderomotive 
drive, and the space overlap of these two fields, matter, we conclude that accounting for a 
diffraction term in the EPW envelope equation is not essential. We actually checked this by 
including a linear diffraction term in the EPW envelope equation and found no significant 
difference in our results for Raman reflectivity. 

Now, if we denote E px = Eoe llfx , with Eq = \E px \, then, provided that the plasma wave 
and the laser drive are perfectly in phase, Eq. (IB.7j) yields 

OE E 

+ v g0 (k p /k )VE + Vg0 ^-V x .k p + u hn E Q = TplE^l, (B.8) 

where k p = k x + V(p x . Eq. ( 1B.8[) is just the linear limit of Eq. ( ITS"]) , except that it is 
written for the x-component of the field amplitude and not for its norm (actually, replacing 
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Eq with E p only yields terms of higher order than that at which the envelope equations were 
derived). We therefore conclude that Eq. (1161) derived in Section [Til only misses the phase 
modulation due to the transverse variation of the field amplitude. More precisely, still when 
the EP W and the laser drive are perfectly in phase, Eq. (IB. T|) also yields 

d k ir [SK + (6k 2 y + 5k 2 z )E /2k - (A ± E )/2k ] + dju = 0, (B.9) 

where 8k = Vip x and 5u = —dttp x , which shows that the central wave number and frequency 
of the EPW are, respectively, k p = k x + Sk and u p = u + 5u such that 

1 + £ r (|4U P )] E = d k £ r (A ± E )/2k . (B.10) 

The former equation makes explicit how the EPW dispersion relation is affected by the 
transverse variations of the EPW amplitude. However, in the nonlinear regime, the trans- 
verse phase modulation of the plasma wave is mainly induced by the adiabatic frequency 



shift, 5u) p . Indeed, as discussed in Ref. 



11| . Sojp may by calculated by making use of the 



adiabatic approximation, and therefore at 0-order in the transverse variations of the EPW, 
provided that v t h/l± < w pe /20, where l± is the typical transverse gradient length of the 
EPW and v t h the electron thermal temperature. For the parameters used in this paper, 
Vth/l± ^$ Wp e /20 if lj_ > A;, where A; is the laser wavelength. 
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